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ABSTRACT 

Some formation scenarios that have been put forward to explain multiple popu¬ 
lations within Globular Clusters (GCs) require that the young massive cluster have 
large reservoirs of cold gas within them, which is necessary to form future gener¬ 
ations of stars. In this paper we use deep observations taken with Atacama Large 
Millimeter/sub-millimeter Array (ALMA) to assess the amount of molecular gas within 
3 young (50 — 200 Myr) massive (^ 10® Mq) clusters in the Antennae galaxies. No 
significant CO(3-2) emission was found associated with any of the three clusters. We 
place upper limits for the molecular gas within these clusters of ~ 1 x 10® Mq (or 
< 9% of the current stellar mass). We briefly review different scenarios that propose 
multiple episodes of star formation and discuss some of their assumptions and impli¬ 
cations. Our results are in tension with the predictions of GC formation scenarios that 
expect large reservoirs of cool gas within young massive clusters at these ages. 

Key words: globular clusters: general - galaxies: star clusters: general - galaxies: 
star clusters: individual: W32187, W32604, W31730 


1 INTRODUCTION 

Since the discovery of the evidence of multiple stellar popu¬ 
lations in the colour-magnitude diagrams (CMD) and chem¬ 
ical abundances of globular clusters (GCs), the classic view 
of GCs as si mple stellar populations has been ca lled into 
question fcf. ICratton. Carretta fc Bragaglial l2012l l. Differ¬ 
ent scenarios have been put forward to account for these 
“anomalous” observed features. 

The majority of the proposed formation scenarios for 
GCs adopt multiple generations of stars in order to explain 
the observed anomalies. Essentially the main concept is that 
the chemically processed ejecta of some kinds of stars from 
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the first generation {polluter stars) breed a second genera¬ 
tion of stars with the observed “anomalies”. Several kinds 
of stars have been suggested to be the po lluters namely: 
Asymptotic Giant Branch (AGB) stars (e.g . iD’Ercole et al.l 
I2OO8I : l^nzini|[200l : IConrov fc Spergelll201ll - hereafter DOS, 
R08 and C Sll respectively), fas t rotating massive stars, 
ERMS (e.g. iDecressin et al.ll2009l : iKrause et al. 201^)^_and 
massive stars in interacting binary systems ( de Mink et al.l 

I2OO9II . 


In all of these scenarios, the ejecta from the polluters 
is expected to cool and sink into the gravitational potential 
well of the cluster where it will spawn the future generations 
of stars. These models all require additional “pristine” ma¬ 
terial (i.e. the gas from which the first generation was born) 
to be (re)accreted by the clusters from their surroundings in 
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order to match the chemical patterns found in GCs, or al¬ 
ternatively that some left over gas from the first generation 
stays bound in the centre of the cluster and mixes with the 
pollut ed material and su bsequently forms stars (e.g. DOS; 
CSll; iKrause et al.ll201^ '). 

Each type of polluter star is expected to have a different 
timescale for the second episode of star formation. For exam¬ 
ple, the AGB scenario predicts that the second star forma¬ 
tion episode begins ~ 30 Myr after the birth of the first gen¬ 
eration (when the super-AGB stars begin to eject their low 
velocity winds). This then continues until ~ 80 — 100 Myr, 
when it is truncated by prompt SNe la explosions. For this 
model this is necessary in order to match the observed abun¬ 
dance trends in GGs (e.g. DOS - as lower mass AGBs stars, 
with longer lifetimes, have different abundance yields). How¬ 
ever, other authors have suggested that this maximum needs 
be pushed to later times (~ 200 — 300 Myr) in order to al¬ 
low the Lyman-Werner flux density to drop sufficiently for 
the gas to cool and form new stars (GSll, see also R08 for 
an alternative long timescale AGB scenario - few 10® yr). 
On the other hand, the FRMS and the interacting binaries 
scenarios work on much shorter timescales (a few Myr) for 
the formation of the second generation. 

In this series of papers, we have used observations of 
you ng massive clusters (YMG s) often seen as young GGs 
fe.g. ISchweizer fc SeitzerlllQQSfl . to probe the different mod¬ 
els for GG formation. Some of these studies were focused 
on the prope r ties o f the gas within YMGs. For example 
iBastian et al.l ll2013ll analysed over 100 YMGs for evidence 
of ionised gas ([Olll] or H/3) due to ongoing star formation, 
and found none. The authors set an upper limit of 1-2% 
of the cluster stellar mass to be currently forming within 
the cluster. In a following work, IBastian fc StradCT (l2014l l 
placed limits on the amount of Hi and/or dust found within 
13 LMC/SMC YMGs (with ages between 15 and 300 Myr, 
and masses between 10"* and 2 x 10® Mq) to be < 1% of 
the stellar mass, this result is at odds wit h the expec t ation s 
of some scenarios (e.g. CSll). In addition. iLongmord (l2015l f 
has shown that high levels of extinction/dust would be ex¬ 
pected for the clusters if there were to be such large reser- 
voirs of gas/dust within these clust e rs as suggested by DOS. 
IBastian. Hollvhead fc Cabrera- ziiil (l2014l l also used YMGs 
to test a prediction from the FRMS scenario, namely that 
massive young clusters should remain embedded in the natal 
gas cloud for the first 20 — 30 Myr of the cluster’s life. How¬ 
ever, it was found that YMGs, independent of their mass, 
are very efficient at expelling any gas within/near them, be¬ 
ginning 1 — 3 Myr after their formation. Similar results have 
been found in YMGs in the Antennae, where they appear to 
have removed the molecular gas within them over the first 
10 Myr of their life clearin g the gas out to a radius of ^200 
pc ilWhitmore et al.|[20l4l . 

In the current paper we look for evidence of molecular 
gas necessary to form the second generation of stars within 
young GCs with Atacama Large Millimeter/sub-millimeter 
Array (ALMA) observations. Our study is centred on 3 
young (< 200 Myr) massive (-^ 10® Mq) clusters in the 
Antennae galaxies (W32187, W32604 and W31730). We es¬ 
timate the escape velocities for these clusters to be between 
50-130 km/s, well above the velocity of the stellar ejecta 
of the different polluters, allowing them in principle to be 
able to retain the ejecta from these stars. These clusters 


have very little reddening, with E{B — V) ranging from 0.06 
to 0.5 mag, and are amongst the 12 brightest/most mas¬ 
sive clusters in the Anten nae galaxies, which ma kes them 
ideal for this kind of study dWhitmore et al]|2010l . hereafter 
WIO). 

The paper is organised as follows: In 321 we present the 
ALMA data of these clusters and in jJSIwe show the proce¬ 
dure used to estimate the H 2 masses for these clusters. We 
discuss our results and present our conclusions in (jJIa'iid 351 
respectively. 


2 DATA 

We made use of “Band 7” (345 GHz) ALMA observations of 
the Antennae galaxies (NGC 4038/39) from Cycle 0 project 
2011.0.00876. A detailed descripti on of these observation s 
and data reduction can be found in lWhitmore et ^ (l2014h . 
along with a formal introduction to the highest spatial res¬ 
olution CO (3-2) observations of the overlap region in the 
Antennae. 

In Fig. m we present a colour composite image of the 
clusters of our analysis (based on HST imaging); W32187, 
W32604 and W31730 (from east to west respectively, W10)Q- 
The clusters of interest are highlighted with green circles. 
The contours indicate the integrated CO emission, from 
which it can already been seen that no CO emission is asso¬ 
ciated with the clusters. 

Only these clusters fulfilled the following selection cri¬ 
teria 1) be located in on our ALMA held of view and away 
of any background CO emission 2) were young (< 200 Myr) 
and massive (> 10® Mq) 3) have near-IR spectroscopy used 
to determine radial velocities and velocity dispersions. There 
are 4 more clusters that satisfy criteria 1) and 2) but they 
lack near-infrared spectroscopy (i.e. with no radial velocities 
nor velocity dispersions), these clusters will be studied in a 
future work. 

Fig. [2] shows the spectra of the (primary beam cor¬ 
rected) CO (3-2) emission towards the three clusters in units 
of mjy/beam. For these observations the channel spacing is 
5 km/s and the synthesised beam size is 0.56” xO.43”. All 
spectra were extracted from single beams taken from the ap¬ 
proximate I-band intensity peak. These positions are given 
in Table [T] together with the characteristics of each cluster 
(see 33 ) 1 . We noted a systematic offset of -1.088” and 1.65” 
in R.A. and Dec. respectively, in the positions of the clusters 
reported in table 7 in WIO with respect to the ones in the 
HST images. 


3 MOLECULAR GAS AND STELLAR MASS 
ESTIMATES 

Following equation (3) in iBolatto. Wolfire fc LerovI (l2013l 'l: 

M = 1 n'^xin'^ r _ ^co _ ScqAvD 

mol 1^2 X 1020 cm-2(K km/s)-i J {1 + z) 

we derived the H 2 gas masses for these clusters. For the 
H 2 mass estimates we assume a CO (l-0)-to-H2 conversion 

^ a.k.a. cluster S2_l, [W99]15 and Sl_2 in iMengel et al.l 1 I 2 OO 8 I I. 
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Figure 2. ALMA spectra of the 3 young massive clusters in the Antennae galaxies. The spectrum of each cluster is presented with the 
highest resolution (bins of 5 km/s) and also binned to the approximate velocity dispersion reported in M08. There is no detection of any 
significant CO (3—2) emission in any of the clusters studied in this work. The vertical dashed lines represents the radial velocities of each 
cluster from M08 corrected to LSRK. The position of the extracted spectra and rms values of the flux are reported in Tabled 
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Figure 1. Zoom to the region containing the clusters. This com¬ 
posite image was made with HST observations where the red, 
green and blue colours are from Pa/3, 7-band, R-band images 
ffrom lwhitmore et alj[2014h . The contours overlaid in yellow are 
CO (3-2), representing flux levels of 4, 24, and 48 K km/s. The 
green circles highlight the clusters of interest. 


factor of A'co = 2 x 10^^ cm“^(K km/s)“^, a luminosity 
distance of Dl = 21.8 Mpc to the Antennae galaxies and a 
redshift of 2 : = 0.005688. Finally Sco^v, is the integrated 
CO (1-0) line flux density, in our case given by 


ScoAv = IhfA! 

rsi 


where -Frms is the rms value of the flux of our ALMA CO 
(3-2) spectra shown in Fig. [2] (we assumed this value to be 
the standard deviation of the flux, a, for these spectra); Aw 
is the channel spacing of the spectra (5 km/s for the high 
resolution spectra and 10 or 20 km/s for the low resolution 
ones, see Fig. [2|; and rai is the CO (3-2) to CO (1-0) line 
ratio. We adopted rsi = 0.792, this value was obtained av¬ 
eraging the mea surements of r^i report e d for the Antennae 
on Table 2 from IZhu. Seaauist fc Kunol (1200311 . 

In Table [T] we present the characteristics of these clus¬ 
ters. The ages and photometric masses are from WIO. We 
remeasured the sizes of these clusters and used these value s 
to recalculate the dynamical masses of iMengel et al.l ll2008ll , 
hereafter M08. For the sizes, we made an empirical PSF 
based on AC S F-band drizzled images and used ISHAPE 
llLarsenlll999ll t o measure the effecti ve radii adopting differ¬ 
ent profiles fsee lBastian et ^l2009l for more details on the 
technique). The derived sizes were 2.7, 1.7, and 2 times larger 
than reported in M08 for W32187, W32604 and W31730 re¬ 
spectively, which corresponds to an increase of the dynami¬ 
cal mass of the same proportion. We find excellent agreement 
between the photometric and dynamical masses. 

In this table we also include Elms, ScoAv and H 2 mass. 
We note that the measurements presented here are Ict upper 
limits for ScoAv and H 2 mass. Already with this data we 
do not find clear evidence of molecular emission above the 
la level in any of these spectra. 
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Table 1. Clusters’ ages, spectra position, photometric masses, projected half-light radii, dynamical masses, flux spectra rms, Scq An and 
H 2 mass. Values in parenthesis correspond to the binned spectra. 


Cluster 

a 

<5 

log(Age)“ Mph“ 

rhp 

'^dyn Frms 

S'oo^'^ H 2 mass 

ID 



(xlO® Mq) 

(pc) 

(xlO® ^ 0 ) (mjy/beam) 

(mjy km/s) (xlO^ ^o) 


W31730 

12 01 55.462 

-18 52 19.108 

7.7 

1.7 

4.1 

1.4 

3.98 (3.21)*’ 

25.11 (40.58)*’ 

1.25 (2.01) 

W32604 

12 01 55.354 

-18 52 17.824 

7.7 

1.6 

3.4 

2.4 

4.48 (2.84)“ 

28.29 (71.73)“ 

1.40 (3.56/ 

W32187 

12 01 55.184 

-18 52 18.442 

8.3 

1.9 

8.1 

1.8 

4.43 (3.47)*’ 

27.98 (43.84)*’ 

1.38 (2.17) 


“ from WIO. We assumed conservative errors for these age estimates of ±50%. 

^ value for the spectrum binned to 10 km/s (similar to the velocity dispersion of 11.5 km/s reported in M08 for W31730 and W32187). 
value for the spectrum binned to 20 km/s (similar to the velocity dispersion of 20.2 km/s reported in M08 for W32604). 


4 DISCUSSION 

From Table [T] we can see that our non-detections of CO 
in the ALMA spectra of these clusters can be translated 
into upper limits on the H 2 mass estimates which are of 
the order of ~ 9% of the stellar mass. Having said that, 
we note that this estimate is subject to the choice of Aco- 
For this stu dy we have used the value reporte d for the Milky 
Way disk bv lBolatto, Wolfire fc Ler^ (l2013ll . However, sev¬ 
eral studies have shown that the Xgo in the Antennae is 
lower by a fa ctor of 2-10 (e.g. IZhu. Seaauist fc Kund[2003l : 
iBisbas et al.ll20r3 l than the Milky Way value. As a conse¬ 
quence the H 2 mass estimate will decrease by the same factor 
if Ago decreases, so the values of moleeular gas reported in 
Table[^are striet upper limits. 

As mentioned in m given the age range of these clus¬ 
ters (50-200 Myr) we are only able to place constraints on 
the AGB scenario. However, this sample spans the critical 
timeframe where the second generation of stars is expected 
to be formed for this choice of polluters (30-300 Myr, see 

There are several variations of this scenario that adopt 
AGB stars as the source of enrichment. In the following sub¬ 
section we outline the essentials of an “AGB scenario” and 
some distinctions between different scenarios. 


4.1 AGB scenario: The essentials 

(i) A first generation (IG) is born as a simple stellar pop¬ 
ulation, i.e. all the stars with the same age, and homoge¬ 
neous (pristine) chemical composition. 

(ii) The remaining pristine gas that formed the IG, is 
expelled from the cluster at some time during the following 
~ 30 Myr by the winds of massive stars and core-collapse 
SNe (R08). 

(hi) A second generation (2G) starts to form. The first 
stars to be born will be from the pure (i.e. undiluted) ejecta 
of the most massive (~ 8 —9 M 0 ) AGB starfl This material 
will be forming the 2G stars with the most extreme abun¬ 
dance patterns, that is, the most He-rich population which 
is also the one with the most 0-poor stars observed today 
in old GGs. We will refer to this p opulation as the extre me 
2G (the extreme stars identified bv ICarretta et aDl2009l j. 
Guriously, the 2G stars are assumed to span a rather narrow 
mass range 0.1 — 0.8 Mq (DOS). This preference for a trun¬ 
cated initial mass function (IMF) for the 2G stars in these 


^ a.k.a. super-AGB stars 


models, would have two main consequences 1) this genera¬ 
tion of stars will not have type H SNe that would clear the 
gas, interrupting the formation of the 2G; and 2) minimise 
the “mass budget problem”, i.e. the fact that one would re¬ 
quire a IG many times more massive than the amount of 
stars with primordial abundances seen today. Both conse¬ 
quences would be treated again in (vi) and (vii) respectively. 

(iv) Only after the extreme 2G is born, the cluster starts 
accreting pristine material from a surrounding reservoir. The 
second star formation episode will continue forming stars, 
this time with intermediate chemical abundances (i.e. abun¬ 
dances between the extreme and pristine ones) from the 
ejecta of the lower mass (< 8 Mq) AGB diluted with this 
accreted material^ 

We need to underline this point, given that the dilution is re¬ 
quired in all the scenarios that share AGB stars as polluters. 
Although several attempts have been presented, to date it 
still remains quite equivocal where was this pristine material 
stored, what triggers the accretion, what process regulates its 
rate, and what mechanism ends it. However, it is possible 
to constrain the nature of this sequence of events with the 
current knowledge of the field, for example: 

• Where does this pristine gas reservoir come from? 

This element is not clear in the scenarios so far. DOS 
implies that the IG cluster is surrounded by some pris¬ 
tine gas, and the accretion starts after the extreme 2G is 
formed. On the other hand GSll, argue that the young 
GC continuously accretes material, as it sweeps the sur¬ 
rounding interstellar medium (ISM) in a competing pro¬ 
cess against ram pressure. Later, after the number of en¬ 
ergetic photons from the IG drops, the gas that has been 
accumulated in the cluster’s potential is allowed to cool 
down t o form stars (CSll). _ 

• In lD’Ercole, D’Antona fc Vesperinil ((2oT3) (hereafter 
Dll), the authors conclude that the maximum amount 
of pristine gas that contributes to the 2G formation is 
at most equal to 10% of the IG mass in order to re¬ 
produce the extent of the observed abundance patterns. 
This limit was obtained assuming high star formation ef- 

® CSll in their model predicts a continuum accretion of pristine 
gas just after the IG is formed, but the star formation is inhib¬ 
ited during the first 100-200 Myr given that the Lyman-Werner 
IG photons have been keeping the gas warm. Due to this (ii) 
and (iii) never took place in this variat ion of the AGB scenario. 
iD’Ercole, D’Antona fc Vesperini ll201l!) argue that this mecha¬ 
nism (continuous accretion) would not be able to reproduce the 
extreme populations observed in GCs, see E31 
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ficiencies (~ 100%), if lower values were in play in the 2G 
formation, the mass of the gas reservoir should increase 
(further discussed in 1)4.2|l . 

• The expulsion of the pristine gas mentioned in step 
(ii) by core-colapse SNe, should have been through a very 
distinct mechanism in the sense that while it is able to 
expel this pristine material (now polluted with Fe-peak 
and a-elements) away from the cluster, it does it in such 
way that this polluted material never reaches the pristine 
gas re servoir. This has been explored by R08 and iRenzinil 
(l2013l i. concluding that if only 0.2% of the core-collapse 
SN ejecta would have reached the pristine gas reservoir the 
2G of the cluster would exhibit [Fe/H] inhomogeneities 
like the ones observed in ojGen, M22 and Terzan 5. 

It has been put forward that some asymmetries in the 
gas distribution could prevent the Fe-polluted material to 
reach the reservoir, avoiding the contamination (cf. DOS). 
In this case the core-collapse SNe leave an “hour glass” 
cavity after clearing the leftover gas from the IG, allowing 
the uncontaminated material around the “hour glass” to 
collapse back into the cluster for step (iv). 

• The accretion must be an accelerated process. The 
theoretical yields of the AGB stellar models predict a cor¬ 
relation between Na and O, in their ejecta, as a func¬ 
tion of stellar mass/time. Given this, after the extreme 
2G is formed, the pristine gas accretion rate should in¬ 
crease rapidly in time, in order to be able to dilute the 
material not just compensating their natural correlation 
trend but to a higher degree, in order to reproduce the 
observed Na-0 anti-correlation. 

(v) It is assumed that the star formation of the 2G took 
place in bursty episodes. This premise in the models is 
backed by the evidence of discrete (i.e. non continuous) He 
abundances in the subpop ulations of some mas si ve GG s (e.g. 
NGG 6397, 2808 & 6752 - iMilone et aLll2012al l^. [20ili . Fur¬ 
thermore, this assumption also agrees with the recent dis- 
covery of a few massive GGs (e.g. NGG 6752 & 2808 cf. 
iGarretta et al.ll2012l : ICarrettall2014h showing clear discrete¬ 
ness on their light element abundance patterns, specifically, 
the Mg-Al anti-correlation. 

(vi) The formation of the 2G should end before the low- 
mass (< 3 Mq) AGB stars produce large amounts of G, 
whic h is especially the ca se for low metallicity AGBs (cf. 
R08. iD’Ercole et al.ll201ol l. The production of C, would af¬ 
fect the sum of GNO elements which appears to be constant 
(within a factor of 1.5) among primo rdial and enriched 
stars in GGs (e.g. IGarretta et al.l [200511 . This requieres a 
mechanism which should be synchronised to start operating 
just before these lower mass AGB go into action, otherwise 
(in principle) the star formation would continue. Regarding 
this, most authors have agreed on the role of prompt (~ 10® 
yr) la SNe on this matter. 

In CSll, the 2G star formation happens suddenly. Here 
an instantaneous burst of star formation gives birth to the 
2G of stars after the number of Lyman-Werner photons 
( 912 A < A < llOOA) drops sufficiently to allow the star 
formation (200-300 Myr after the birth of the IG, GSll). 
The star formation ends when the type II SNe from the 2G 
and prompt la SNe from the IG start to go off, keeping the 
cluster gas free for the rest of its life. In a similar manner, 
for R08 the formation of the 2G is likely span up to 300 


Myr as well. Gonstrastingly, in D08 the prompt la SNe are 
expected to end the 2G formation episodes ~ 80 Myr after 
the IG star formation event. In this model there will be no 
type II SNe from the 2G IMF since all the stars have masses 
well bellow the limit for core-collapse SNe (D08). 

(vii) After the 2G is formed in these models, the young 
cluster ends up with just a small fraction of the stars 
showing the chemical signature of the AGB material. In 
order to match the observed high fraction (close to 1:1) 
of enriched/non-enriched stars (i.e. ratio of “first-to-second 
generations of stars”) these scenarios assume that GGs un¬ 
derwent strong mass-loss and lost most of their first genera¬ 
tion stars, usually 90% or more, which means that GGs were 
singificantly more massive at birth than seen today. Further¬ 
more, these models usually assume that the 2G stars had a 
truncated mass function. 

In D08, the authors claim that the cluster should have been 
about 10 times more massive at birth, however, this value 
has some severe underlying assumptions. One can see from 
Fig. 4 of D08, that the stellar mass of the 2G is just ~ lO"*'® 
Mq, this corresponds to a ~ l/30th of their IG mass (10® 
Mq). So only if the current (old) cluster has a particularly 
low fraction of enriched/non-enriched stars (1:3) and only 
if IG stars were lost (all 2G stars were kept in the cluster) 
one can reach this value for the mass of the IG (originally 
10 times more massive). Alternatively, if the current clus¬ 
ter has an even fraction (1:1) of enriched/non-enriched stars 
the cluster should have been 30 times more massive at birth 
(again, this value assumes no 2G stars escaped the cluster 
during its dynamical evolution). 

Additionally, as mentioned in (iii), D08 is assumed a trun¬ 
cated IMF (only forming stars with 0.1 — 0.8 Mq) for the 
2G. If the 2G would have formed following a conventional 
(i.e. lKroupa. Tout fc Gilmorelll99^ non-truncated IMF (i.e. 
stars form with masses between 0.1-100 Mq) the mass of the 
first generation should scale by a factor of 2 (i.e. IG should 
have been 20 or 60 times more massive in the past depend¬ 
ing of the assumed fraction of enriched/non-enriched stars 
of 1:3 and 1:1 respectively). 

On the other hand, in GSll the young GGs had a similar 
masse s than the ones observed today (Dll), while l^onrovl 
(I 2 OI 2 II support that the GGs were at least 10-20 times more 
massive at birth. 

Although there has not been reported any evidence of 
a strongly non-standard IM F in GGs or YMCs (cf. 
iBastian. Covev fc Meveill201Cll ~) this option is difficult to ex¬ 
clude. However, in the last few years there have been some 
constraints reported regarding the strong mass-loss. 

_ Larsen. Strader Sz Brodi3 (l2012l ~) and iLarsen et al.l 

1 I 2 OI 4 I 1 studied three dwarf galaxies, and looked at the 
amount of GGs stars which could have been lost and now 
form part of the field populations of their host galaxies. In 
both studies the authors found high GCs-to-field stars ratios, 
concluding that these GGs could only have been at most 4-5 
times more massive initially. These values are upper limits, 
given that it was assumed that all of the field stars that 
share the same metallicity as the GGs (in a broad metallic¬ 
ity range), formed in these clusters (i.e. no field stars formed 
with similar abundances) nor did they take into account that 
lower mass clusters would likely have be en formed at th e 
same time and have since disrupted (e.g. lKruiissenll2015l L 
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contributing t o the field popu lation. Furthermore, the re¬ 
cent model bv iKruiissenI (l2015l i. putting GC formation into 
a hierarchical galaxy assembly context, also found that clus¬ 
ters that survive to the present epoch were likely to be only 
2-3 times more massive than currently observed. 

In the following subsection we will present a “toy 
model” designed to explore the viability and consequences 
of the AGB scenario, given the constraints we have just pre¬ 
sented regarding the stellar mass and upper limits on any 
potential gas within the observed clusters. 


4.2 Toy model 


The “toy model” is based on the most massive young globu¬ 
lar cluster in our sample, W32187 with a mass of ~ 2 x 10® 
Mq. This cluster lies within the age range where a second 
episode of star formation is expected to occur in the AGB 
model. We adopt the initial mass up p er limits discussed 
above jLarsen. Strader"^ Brodid l20ld : iLarsen et al.l l2014l : 
iKruiissea 2015l l. and assume that the current mass of the 
cluster is 4 times more massive than it will be in ~ 10 
Gyfl i.e. assuming that the current cluster stellar mass is 
the initial mass. This means that as an old GC, this cluster 
will have a mass of M* = 5 x 10® Mq. If we also assume a 
1:1 ratio for the number of £rst-to-second generation stars at 
an age of ~ 10 Gyr, this would correspond to a mass of the 
second generation (2G) stars of MJ® = 2.5 x 10® Mq. This 
is a lower limit to the amount of gas that would need to be 
present, integrated over time. However, if we adopt a stan¬ 
dard starTormu^nefficiency for this cluster of SFE=l/3 
fe.g. iLada fc La(5l2003h . this means that Mg^^ = 7.5 x 10® 
Mq of molecular gas is required to give birth to the 2G. 
The rest of this gas i.e. the other 2/3 that will not end up 
in stars, is expelled from the cluster (see 1)4.511 after the star 
formation takes place. Given these assumptions, how much 
molecular gas would we expect at any given time within this 
cluster? 

To answer this question we consider the window where 
star formation can take place according to DOS, At = 50 Myr 
(i.e. from 30-80 Myr, see CT . If the star formation is ex¬ 
tremely rapid, in the sense that any gaseous material within 
the cluster is used in the star formation very quickly, e.g. 
At|® = 1 Myi0, we would expect to find Mg®s = 1.5 x 10^ 
Mq (i.e., 1/50*® of the total gas mass) at any given mo¬ 
ment during this 50 Myr time window, which lies below our 
detections limits of ~ 1 x 10® Mq. However, if the star for¬ 
mation event lasts longer, i.e. that any gas within the cluster 
is present for say At^ — 10 Myr, before being consumed in 
star-formation, the amount of gas at any given time would 
increase to Mg)((, = 1.5 x 10® Mq, which is above our (con¬ 
servative/overestimated) detection limit for this cluster of 
1.38 x 10® Mq. Moreover, if such clusters have undergone 
a significant mass loss (as claimed in the D08 model), then 
we would not be looking at the initial mass, but rather the 


^ This value differs the original DOS scenario (as seen in 
mi. however, with this adjus t ment we take into a c count the 
Ibarscn. Strader fc Brodi3 ll2012ll : iLarsen et al.l ll2014h : iKruiiss^ 
ll2015l ~) results. 

® We note that if we assume a slighter shorter timescale of 
AtgJ, = 0.6 Myr our simple model is in good agreement with 
DOS predictions. 


present day values (contrary to what was assumed here). 
Hence, if the initial mass was significantly higher, we would 
expect more gas to be present in the cluster, for example 
if the cluster has already lost 50% of its initial mass, then 
these estimates would increase by a factor of 2. 

CSll suggest that must be long (~ 100—300 Myr) 
as the Lyman-Werner flux density from the first generation 
stars will be too high to allow the gas/dust to cool suffi¬ 
ciently to form stars. In this case, we would have expected 
to detect large amounts of gas (> 10% of the stellar mass) 
within all three of the clusters observed in the present work, 
in contrast to what is observed (an upper limit of < 9% of 
the stellar mass). 

This simple “toy model” demonstrates the general 
bounds that may be expected in the AGB scenario. The 
actual amount of gas within the cluster depends critically 
on 1) the assumed SFE and 2) the time of which the ISM 
may exist within the cluster before being expelled or used 
to form 2G stars. 


4-8.1 Summary of the toy model and comparison with 
previous works 

Our “toy model” was developed to find the minimum 
amount of gas during the formation of a 2G needed to form 
the number of stars with enriched abundances (a.k.a. 2G 
stars in the AGB scenario) seen today in a GC. In this model 
we have assumed that at every Atg®, 1/3 of all gas (regard¬ 
less of its origin i.e., either ejected by AGB stars or accreted 
from the pristine gas reservoir) is turned into stars, while the 
remaining 2/3 of the gas will be lost from the cluster (this 
could be by means of any of the mechanisms mentioned in 
1)4.511 . Nevertheless, this is just an assumption. However, if 
this material is not cleared away from the cluster, but in¬ 
stead remains within the cluster, it could build up over time. 

This alternative has been modelled and thoroughly de¬ 
veloped by D08. Here, the authors assume that the gas that 
has not found its way into a star (i.e. 2/3 of the total gas 
in our case) remains in the cluster. Perhaps this material is 
too warm to form stars yet, but the key issue here is that 
it has not been expelled from the cluster (in D08 this gas 
will eventually cool and start forming new stars j3- This pro¬ 
cess results in the accumulation of gas within the centre of 
the cluster. However, since the gas mass and density are in¬ 
creasing, DOS assume that this will also have an effect on the 
star formation rate (SFR), with higher densities resulting in 
higher SFRs. This treatment is physically motivated, but of 
course adds some additional assumptions and parameterisa- 
tions that are beyond the scope of our simple toy model. 

However, if some kind of feedback is considered (which 
is not the case in DOS - see >)4.51l and would have kept the 
gas from cooling (so that it cannot form new stars, but it 
does remain in the cluster), then the gas mass will increase 
as a function of time inside the cluster (e.g., as in the CSll 
scenario). The longer the gas accumulates within the cluster 
the easier it would be to detect it. Additionally, if the gas 


® If the gas can cool, and remains in the cluster, it will eventually 
form new stars. Hence, even if the SFE is low per unit time in 
this model (DOS), it can approach unity over the full formation 
epoch of the 2G (i.e. assumed to be 50 Myr in this case). 
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temperature increases (due to the feedback mechanism, e.g., 
Lyman-Werner flux) then the gas would be easier to detect, 
for a given amount of gas mass, with the type of observations 
used in the current work. Although perhaps our toy model 
might be considered less physically motivated than assuming 
some fraction of gas is retained, we note that this approach 
provides a lower limit to the amount of gas expected in the 
cluster at any given time. 

Finally, we must remark that our unsophisticated model 
is i n no way meant to replace the ones provided by 
DOS/ lD’Ercole et al] ll2010ll or CSll which are much more 
physically motivated, which is the reason we why we have 
chosen to compare against them. Nevertheless, we have 
shown that this simple description in our “toy model” is 
able to reproduce the same processes as the more elaborate 
model of DOS, without the need to include further uncon¬ 
strained parameters. Here we have shown that a good agree¬ 
ment is achieved between DOS model and our “toy model”, 
if one adopts a relatively short At^ of 0.6 Myr. 

4.3 Accretion of pristine gas 

CSll have presented a model for how pristine gas may be 
continuously accreted at a (nearly) constant rate by the GC 
in order to dilute the AGB ejecta within the cluster, and 
form 2G stars ( such strong dilution is required on chemi¬ 
cal grounds, e.g. iD’Ercole et al.ir2010h . In the CSll model, 
young GCs contain 10% of their initial stellar mass in gas 
within the cluster, and this gas acts as a net as the cluster 
passes through the galactic ISM, picking up further pristine 
material. In this case, the clusters would always be expected 
to have > 10% of their stellar mass in gas within them, which 
would be above our detection limit. Hence, this scenario is 
not consistent with the observations presented here. 

Moreover. ID’Ercole. D’Antona fc Vesperinil (I 2 OIIII have 
shown that any model with nearly constant accretion of pris¬ 
tine gas can not reproduce the observed abundance trends 
found in GCs. Instead, the accretion of pristine gas must 
be finely tuned, i.e. it cannot contribute significantly until 
~ 20 Myr after the 2G begins to form. This allows for the 
most enriched 2G stars to form entirely from (super)AGB 
ejecta, while subsequent 2G stars form from a combination 
of AGB ejecta and pristine gas, which dilutes the abun¬ 
dances, and turns the predicted correlation between Na-0 
from AGB stars, into an anti-correlation (as mentioned in 

EH) 

4.4 DOS mass prediction 

From Fig. 3 of DOS, we obtained an estimate of the amount 
of gas expected to be found in a young (100 Myr) massive 
(10® Mq) cluster, integrating the gas density profile in a 
radius of 25 pc (approximately the size of our beam read- 
ius). We obtained ~ 1500 Mq of gas, which corresponds to 
a reservoir of molecular gas of ~ 0.15% of the stellar mass 
at an age of 100 Myr. This implies an extremely rapid 
(Atgf ~ 6 X 10®yr, see 114.211 . As discussed in CSll, such 
rapid turn-over from accreted/ejected gas to stars, is likely 
optimistic, as it ignores the large heating effect of evolving 
first generation stars, i.e. stellar feedback. 

In DOS the only heating source in the standard model was 


the thermalisation of the kinetic energy of the stellar winds. 
DOS ran additional models with “extra energy sources” (e.g. 
X-ray binaries and planetary nebulae nuclei) and found that 
such heating may terminate any 2G star formation (above 
certain threshold). However, detailed models are required to 
understand the exact efficiency of stellar feedback on any 
material residing in the cluster . Such models have ye t to be 
done for young GCs, however iMcDonald fc Ziilstral (l2015li 
have run these models for older GCs and have found that 
stellar feedback is very efficient at clearing out and gas/dust 
within the cluster. Since the feedback at young ages (< 200 
Myr) is significantly higher than at ~ 10 Gyr, it is likely that 
it can provide an effective gas/dust clearing mechanism, po¬ 
tentially explaining the results presented here. 

Nonetheless, due to the lack of heating included in the cur¬ 
rent DOS model, the accreted/ejected gas can sink in the 
very centre of the cluster (< 0.1 pc), where the gas density 
becomes extremely high, leading to rapid and efficient star- 
formation. Additionally, this model assumes that the 2G of 
stars only make up 3% of the initial mass of the first gen¬ 
eration, which needs to be increased in order to match the 
relative numbers of 1G/2G stars in some clusters. 

4.5 Why are YMCs gas free? 

As mentioned in the previous section, DOS studied the pos¬ 
sibility of models with extra energy sources (only the wind 
thermalisation energy was included in the standard model). 
They concluded that if a critical value of Qcr = 2 x 10®® 
erg/s was exceeded, a 10® Mq cluster would abruptly tran¬ 
sition from a model hosting a cooling flow, which would be 
collecting the AGB ejecta in the centre of the cluster, to a 
model hosting a wind (i.e. expelling the material from the 
cluster). Furthermore, they arrive at the conclusion that if 
a luminosity of Q = 1.38 x 10®^ erg/s would to be reached 
within a massive (lO’^ Mq) cluster, the immediate onset of 
a fast wind would preclude the formation of any 2G0. These 
authors proposed X-ray binaries as a possible source of this 
energy. 

I Power et al.l ll2009h studied the role of high mass X-ray bi¬ 
naries (hereafter, HMXBs) in the first few hundred Myr of 
GCs. In their Fig. 2, one can see the evolution of the H 
ionizing photon rate over time, produced by the HMXB 
population of a 10® Mq cluster. In this plot one can see 
that even if only 50% of the high mass binaries that survive 
the SNe of the primary star of the binary system become 
HMXBs (it is assumed that only 1/3 of binaries will survive 
the primary SNe, i.e. ~ 15% of the initial high mass binary 
population), one would expect to have Q ~ 1.5 x lO'^® erg/s 
and Q ~ 2.2 X 10®® erg/s at 30 and 80 Myr respectively. 
From this we see that, during the timescale where the 2G 
is expected to be form in DOS, these Q values are several 
order of magnitude above the threshold {Qcr ~ 10®® erg/s) 
necessary to stop the star formation of the 2G in these 10® 
Mq clusters. However, a crucial parameter that requires fur¬ 
ther modelling is what fraction of this energy is transferred 
to the ISM within the cluster, and what fraction escapes the 

^ The equivalent value for a 10® Mq cluster was not provided in 
DOS, but we assume it would be below this Q = 1.38 X 10®^ erg/s 
limit. 
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cluster directly. 

This makes HMXBs a potential candidate to explain why 
these clusters appear to be gas free. This possibility is reas¬ 
sured when considering that in these environments a sin¬ 
gle HMXB source has a typical luminosity in excess of 
Lx ^ 10^^ erg/s (e.g. lFabbiano, Zezas fc M urray! I 2 OOII : 


IWolter fc Trinchierill20M: Rangelov et al.r2012l'l. 

Other plausible mechanism able to remove the gas from the 
cluster could be ram pressure striping by the ISM around 
the cluster, as shown in CSll. Originally, this mechanism 
would not be effective for the DOS scenario as it is, due to 
the fact that in this model most of the gas gathered by the 
cooling flow would be concentrated in the innermost region 
(< 0.1 pc) of the cluster. Such an extremely concentrated 
volume of gas would not be sensitive to the influence of 
ram pressure, since this mechanism is heavily dependant on 
the surface area/cross section of the gas within the cluster. 
Hence this effect can be neglected for the current version of 
this model. Notwithstanding, this changes if the gas were to 
occupy a larger volume within the cluster. For example, if 
this gas were to be exposed by some heating (some sort of 
stellar feedback, currently not considered in this model e.g., 
Lyman-Werner photons) it would be reasonable to expect 
the density of this gas to decrease, and as a consequence, 
covering a larger volume within the young cluster. There¬ 
fore, if this (heated)gas would now cover the central 1 pc of 
the cluster, the new surface of this gas would have increased 
by a factor of 100, and the effect of the ram pressure would 
increase as well in the same proportion. To summarise, if 
the gas ejected by the polluters stars were not to be concen¬ 
trated in a very small volume at the centre of the cluster, 
it is likely that the ram pressure striping by the ISM sur¬ 
rounding the cluster would play a non-negligible role in the 
expulsion of some the gas out of the cluster. 

For the moment, we can not provide a definitive answer to 
why are YMCs gas free, and in order to reach a definite 
conclusion on what could be the possible explanation of this 
phenomenon, we require further studies that are beyond the 
scope of the current paper. Nevertheless, we have shown that 
at least these two possibilities could play a role in its answer 
and should be contemplated in future studies. 


5 CONCLUSIONS 

We use ALMA CO (3-2) observations to search for cool gas 
within three young massive clusters in the Antennae merging 
galaxies. Some scenarios for the formation of the observed 
multiple populations in GCs predict that such YMCs should 
contain significant amounts of cool gas (^ 10% of the initial 
stellar mass of the system) within them in order to form 
multiple generations of stars (e.g. DOS, CSll). The observa¬ 
tions of clusters W31730, W32604 and W32187 used in the 
current work provide upper limits of ~ 1 x 10® Mq of H2 gas 
within the three clusters or 7.4, 8.8, and 7.3 % of their stel¬ 
lar mass, respectively. We note, however, that the adopted 
Yco may be overestimated by up to an order of magnitude, 
which would decrease the measured upper limits by the same 
amount. 

We have estimated the amount of gas that is expected 
to be present within the clusters, in various versions of the 
AGB scenario, and found that the exact amount depends 


crucially on a number of assumptions and uncertainties (e.g. 
the amount of first generation stars lost, star-formation ef¬ 
ficiencies, stellar IMF variations), hence we cannot reach 
definitive conclusions regarding the AGB model in general. 


For the D08 scenario, the results from W31730 and 
W32604 clusters are consistent with their prediction. How¬ 
ever, if one were to adopt less extreme values for the main 
parameters in play in the model (e.g. IMF, SFE, fraction 
of enriched/non-enriched stars, IG and 2G mass loss, etc.) 
instead of the adopted by the authors, one would expect 
significant amounts of gas in these young clusters (as shown 
in and would be easily detectable by our observations. 
On the other hand, cluster W32187 (200 Myr) in principle 
would lie outside the timescale for the formation of the 2G 
adopted in this scenario (30-80 Myr) established by the first 
SNe la explosions. Having said that, we should emphasise 
that the timescale of these prompt la SNe is highly uncon¬ 
strained, and in this model there is nothing inherent in the 
80-100 Myr values adopted. Gurrent evidence suggests, that 
these events could start as soon as just ~ 40 — 50 Myr or as 
late as 400 Myr after the IG was formed fe.g. iR uiter e t al.l 
l2013l : iMaoz, Sharon fc Gal-YamI I 2 OI 0 I : iBrandt et al. 2010l L 
Hence, the results from this particular cluster have a con¬ 
straining affect in D08 as well. Ultimately, the GSll AGB 
scenario, predicts gas masses well in excess of the upper lim¬ 
its provided here (during these ages), hence our results are 
in strong tension with this model. 


The lack of evidence of any cold gas within these clus¬ 
ters could be explained in two ways 1) The actual levels are 
below our detection limits or 2) that the clusters are very ef¬ 
ficient at removing the gas. This could be due to heating by 
stellar feedb ack and/or subsequent stripping i.e. ram pres - 
sure IGSll. lPower et al.ll2009l : [McDonald fc Ziilstral[20151 1. 
However, further full radiative transfer calculations are re¬ 
quired to assess the precise magnitude of these effects. 


Although the constraints placed by our ALMA obser¬ 
vations are not conclusive in ruling out the AGB scenario, 
we note that they are in agreement with a growing num¬ 
ber of studies that have not found evidence of extended 
(or multiple) star-formation episodes within massive clus¬ 
ters. For example, no evidence of ongoing star formation 
has been detected (down to limits of few percent of the 
mass of the first generation) in young (10® — 10^ yr ) mas- 
sive (10** — 10^ clust ers (cf. iBastian et al.ll201^ . Also, 
iGabrera-Ziri et ^ ll2014l ') used an integrated optical spec¬ 
trum of a massive (> 10^ Mq) young (~ 100 Myr) cluster, 
in order to estimate its star-formation history. The cluster 
was found to be consistent with single stellar generation, 
with no evidence of a secondary burst down to mass ra- 
tios of 10-20% of th e curr ent cl uster mass. Addit i onally , 
IBastian fc Silva-Villal 1 I 2 OI 3 II and iNiederhofer et all (1201'^ 
have searched for evidence of age spreads in resolved mas¬ 
sive LMG/SMG clusters, and found that all are consistent 
with a single burst of star-formation. Altogether, the evi¬ 
dence from observations of YMGs places strong constraints 
on the AGB scenario, ruling out large areas of parameter 
space. 
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